Dopant-Induced Local Pairing Inhomogeneity in Bi2Sr2CaCu208+5 
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A new theoretical model is presented to study the nanoscale electronic inhomogeneity in high- 
Tc cuprates. In this model, we argue that the randomly distributed out-of-plane interstitial oxygen 
dopants induces locally the off-diagonal (i.e., hopping integral) disorder. This disorder modulates the 
superexchange interaction resulting from a large- [/ Hubbard model, which in turns changes the local 
pairing interaction. The microscopic self-consistent calculations shows that the large gap regions 
are registered to the locations of dopants. Large gap regions exhibit small and broader coherence 
peaks. These results are qualitatively consistent with recent STM observations on optimally doped 
Bi2 Sr2 CaCu2 Og+s- 

PACS numbers: 74.25.Jb, 74.72.-h, 74.50.-|~r, 74.20. Fg 



One of the most striking features of high- Tc cuprate, 
Bi2Sr2CaCu208+5 (BSCCO), is the nanoscale electronic 
inhomogeneity as observed, for example, byt he scanning 
timnehng microscopy (STM) 0, S H li B 0- Earlier 
STM experiments shows that the spectra gap varies dra- 
matically on a very short length scale, typically in the 
scale of 50A. At the same time the low energy part of 
the tunneling spectra are extremely homogeneous spa- 
tially, except in the immediate vicinity of known defects 
in the Cu02 plane Coherence peaks in regions 

with large gap are broadened and their height is reduced 
or peaks are entirely absent. Several scenarios have been 
proposed to understanding this surprising electronic in- 
homogeneity. It is commonly accepted that the hole dop- 
ing in the Cu02 planes is proportional to the concentra- 
tion of out-of-plane oxygen dopant atoms. Thus the out- 
of-plane doping also introduces disorder. Therefore, it is 
tempting to speculate that poorly screened electrostratic 
potentials of the dopant atoms cause a variation in the 
local d oping co ncentration, giving rise to the gap modu- 
lations [liJ, |ll|, . Alternatively, one would ascribe the 
distinct electronic structures to the competition between 
two different electronic orders that exist intrisically in 
the system — In the small gap regions, with high intensity 
sharp coherence peaks, the superconductivity wins; in the 
large gap regions, a pseudo-gap phase dominates with a 
competing order such as spin Q> ^ or orbital 0| 
antiferromagnetism, or charge density wave |l7j . The 
competing order scenario is supported by the STM data 
on the underdoped samples llSl Il9l or on the optimally 
doped ones in the mixed state |2(]j , where a checkerboard 
pattern in the local density of states (LDOS) has been 
observed. One distinctive feature of all these scenarios is 
that a uniform pairing interaction throughout the system 
has always been assumed. 

Recent STM data on the optimally doped BSCCO 
showed that the locations of oxygen dopants are posi- 
tively correlated with the large gap regions with small 
charge density variations. It was proposed by Nunner et 
al. |2l| that the dopant atoms modulate the local pairing 



interaction, i.e., the local attractive coupling g between 
electrons are spatially dependent. In the phenomenologi- 
cal model of |2]J , the modulation of g is suggested to arise 
from a modification of local electronic structure param- 
eters, but no concrete microscopic model was explored. 

In this Letter we propose a microscopic model for the 
pair potential disorder. We argue that it originates from 
the modulation of the local pairing interaction induced 
by the doped oxygen atoms. The specific mechanism is 
based on an extended t-J model with the local superex- 
change interaction being modulated by the off-diagonal 
disorder (i.e., the electronic hopping integral). This spe- 
cific mechanism on how the gap inhomogeneity might be 
triggered by the disorder in p air potential, is also con- 
sistent with the findings |22l l2,l | that ri-type impuri- 
ties rather than ra-type impurities in the Nambu space 
are more relevant in the interpretation of the q struc- 
ture observed in most Fourier-transformed STM experi- 
ments |24^^. 

We find that the off-diagonal disorder, triggered by 
dopant atoms, produces a modulations of the gap across 
the sample. To be consistent with the experiment, we will 
assume that dopant atoms are interstitial. Each dopant 
atom therefore leads to a local modulation of the hop- 
ping matrix element on the neighboring links in Cu-0 
plane. This increased hopping will elevate the local su- 
perexchange interaction, which in turn leads to the local 
increase of the pair potential. Typical length scale of the 
modulation in response to the electronic distortion is set 
by the superconducting coherence length ^ ^ 5a and is 
short in our model. We also find that the gap exhib- 
ited in the LDOS spectra reflects directly the local pair 
potential; and that the LDOS spectra in the regions of 
the large gap have coherence peaks that are moving out 
and are less sharp than the coherence peaks at small-gap 
regions. We believe the effect found here is a generic 
property of any model, where the pairing scale is set by 
the local superexchange. 

As early as in 1987, Anderson suggested the relevance 
of the large-JJ limit of the Hubbard model to the prob- 
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lem of high- Tc superconductivity in cuprates |23 and the 
strong correlations are a central part of the high-Tc prob- 
lem. To attack the present problem, we therefore start 
with a large-?/ Hubbard model with an off-diagonal dis- 
order, defined on a two-dimensional square lattice: 



(ij>- 



Here the quantities c[P annihilates (creates) an electron 
on the ith site with spin a, Uj is the hopping integral be- 
tween between sites i and j, ei is the on-site energy, /i is 
the chemical potential, nia is the spin-polarized number 
operator, and indicates the summation over neigh- 

boring sites. Following the similar treatment to a large-[/ 
Hubbard model for a clean system Is^l , one can 

show that, in the large-?/ limit, the disordered Hubbard 
model is equivalent to the following model Hamiltonian: 
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acting on the subspace of empty and singly-occupied sites 
only, where Si is the Heisenberg spin-i operator, the 
local supcrexchange interaction Jij = Atf^/U. In tradi- 
tional disordered Hubbard model, the disorder is mod- 
eled through a random distribution of on-site (diagonal) 
energy while the hopping integrals are taken to be uni- 
form. This approach cannot generate a direct modulation 
of the local pairing interaction, and the resultant pair po- 
tential can only be modulated near the known impurities 
through a local shift in effective chemical potential within 
a self-consistent solution. 

Here we argue instead that the oxygen dopants mod- 
ulate the local electronic hopping integral, and as a di- 
rect consequence, modulate the local pairing interaction 
through superexchange Jij. For the purpose of this dis- 
cussion, we consider only the optimally doped regime. 
In this regime (and hopefully also in the overdoped 
regime), we believe that the slave-boson mean field ap- 
proach might be a reasonable approximation |3]| . In 
this approach, one first writes the electron operator as 
Cia = blficr, with /io- and 6i being the operators for a 
spinon (a neutral spin-^ fermion) and a holon (a spin- 
less charged boson). Due to the holon Bose condensation 
at low temperatures, the holon operator bi is treated as 
a c-number and the quasiparticles are determined by the 
spinon degree of freedom only. Within the mean-field ap- 
proximation, the Bogpliubov-de Gennes (BdG) equations 
can be derived as 1321: 
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vyf*[l - /(£;„)]} and 5i - 
_ f{En)]} where ks is 



with Hij = -[tij&i&j* -I- 4^Xij] + (ei - Here uV^ and 
w" are the Bogoliubov amplitudes corresponding to the 
eigenvalue En- The resonant-valence-bond (RVB) pair- 
ing order parameters (OP) Aij, the bond OP Xij> and the 
spatially dependent hole density Si are determined self- 
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xa = E„{"r*"j7(i?n) 

l-2Ei,JKP/(i?n)- 

the Boltzmann constant; f{E) = [ex-p{E/kBT) + 1]^^ is 
the Fermi distribution function. The holon c-number is 
determined through bi = ^/Sl. Here without the loss of 
generality, we have assumed bi to be real. With exact 
diagonalization, we solve the BdG equations fully self- 
consistently. Once a converged solution is reached, the 
LDOS at zero temperature is then evaluated according 
to: pi{E) = 2J2n[K\^SiEn-E) + \v^\^S{En + E)] where 
a factor 2 arises from the spin sum. The LDOS pi{E) 
is proportional to the local differential tunneling conduc- 
tance which can be measured in an STM experiment |33l | . 
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FIG. 1: (a) d-wave OP; (b) hole density distribution; (c)-(d) 
X- and y -oriented RVB bond OP. They are obtained from a 
solution to an untied single dopant. 

In the numerical calculation, we construct a superlat- 
tice with the square lattice x Ny as a unit super- 
cell. As detailed in Ref. [s^l, this method can provide 
the required energy resolution. Throughout the work, 
we take the size of the unit supercell Na = 32 x 32, 
the number of supercells Nc = 6x6, the tempera- 
ture r = 0, the Hubbard interaction U ~ 10. As we 
have already mentioned before, the hopping integral 
should be modified by the presence of out-of-plane oxygen 
dopants. Without a detailed first-principle calculation, 
the extent of this modification is unknown. Therefore, 
we parameterize the hopping integral in the form |35l |: 
tij = ^ij + Stij, where t?j is the bare hopping integral. 
We choose to be t = 1 for the nearest-neighboring 
hopping and t' = —0.3 for the next nearest-neighboring 
hopping. For the resultant superexchange interactions, 
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FIG. 2; LDOS spectrum at the "on-dopant" site with a gap 
Ai ~ 0.6 (red line) and in a location far away from the dopant 
with a gap Ao ~ 0.5 (blue line). 





I' 



FIG. 3: LDOS maps Pi{E) in real space for energy at i5 = Ai 
(a) and E = Aq (b). 



only those for the nearest-neighbor pairs are considered, 
with a new energy scale J = At^ /U = 0.4. By assum- 
ing that dopant atoms are interstitial, as suggested by 
the STM observations 0, |2^, we restrict for simplic- 
ity the change in the hopping integral to the four bonds 
forming a plaquette in the Cu-0 plane, at the location 
sharply below each oxygen dopant and take = 0.2. 
No site-diagonal disorder are introduced, i.e., Ci = 0. The 
chemical potential is tuned such that the averaged hole 
density from the self-consistency solution is about 0.14. 
The obtained spatial variation of the d-wave is defined as 
Ad(i) = i[A4(i)-hA-i(i)-Ay(i)-A_y(i)]. Throughout 
the following discussion, the energy is measured in units 
of J. 

a single dopant atom. For the purpose of elucidating 
the proposed picture, this calculation is performed with- 
out anchoring the dopant concentration to the averaged 
hole density. We find that both the local pair poten- 
tial (i.e., d-wave OP, see Fig. Q^a)) is enhanced on the 
sites belonging to the plaquette directly below the dopant 
atom (for brevity, we will call these sites as "on-dopant" 
sites). This modulation occurs locally and the "healing 
length" to recover the bulk value is about the bulk super- 
conducting coherence length ^ ba. Associated with the 



local hopping matrix element modulation, the hole den- 
sity also changes. Directly on the off-diagonal impurity, 
the hole density is depressed and cxihibits the Friedel- 
like oscillations away from the dopant site. It means the 
corresponding Bose-Einstein condensed holon field is de- 
pressed at the dopant site. In addition, Fig.n](c-d) show 
the corresponding changes in the x- and y- oriented RVB 
bond OP, which are also recovering a bulk mean field 
solution at a few lattice spacings. 

The LDOS is shown in Fig. El The LDOS calculated 
at the "on-dopant" sites clearly shows an enhanced gap 
with a l/-shaped bottom. Interestingly, it is found that 
the lost spectral weight at low energy is transferred to 
the band edge rather than to enhance the height of the 
coherent peaks at the gap edge. The "on-dopant" gap 
Ai is larger than the bulk value of the gap Aq. Gap en- 
hancement is about 20% for our choice of parameters. 
Enhancement clearly depends on the extent to which 
we increase a local hopping matrix element. However, 
the present result has been sufficient to convey that the 
gap modulation as measured by the STM experiment is 
proportional to the local pair potential modulation due 
to the change in local hopping integral of electrons. It 
is qualitatively different from the case of a single site- 
diagonal unitary impurity, where the impurity-site LDOS 
docs not show the shift of the gap edge though the d-wave 
OP has been strongly suppressed The contrast be- 
tween locally large gap Ai and the sample bulk gap Aq 
is shown in Fig.|3| For the voltage bias at i? = Ai, the 
only bright spot is on the location corresponding to the 
dopant atom. For the bias at i? = Aq, the bulk of the 
sample is bright and the only dark region is on the dopant 
site. 
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FIG. 4: Same as Fig.^but for a concentration of 14% dopants. 
Here we have assumed that each dopant atom produces one 
hole. The distribution of dopant atoms is random. 

many dopant atoms. To strengthen our proposal of 
dopant-induced pairing interaction modulation, we turn 
to the case of multiple dopants, which is more relevant 
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FIG. 5: LDOS pi{E) map at S = Ai (a) and E = Ao (b) for 
a concentration of 14% dopant atoms. 



to the experimental situation. Here we assume that each 
oxygen dopant contribute one hole to the Cu-0 plane. 
Therefore, we tie the concentration of dopant to the av- 
erage hole density in the Cu-0 plane, i.e., about 14%. As 
shown Fig. 21 the large pair potential regions are locally 
registered to the dopant positions, which are randomly 
distributed. The corresponding hole density map shows 
the variations that are anti-correlated with the large gap 
regions. Local x- and y- oriented RVB OP also are regis- 
tered to the dopant locations. Figure |S1 shows the LDOS 
map at the voltage bias at i? = Ai and Aq, respectively. 
It demonstrates again that the intensity map at i? = Ai 
exhibits a positive correlation with the locally large pair 
potential regions (the mixed color areas in Ela)) while 
that at £■ = Ao reflects the small pair potential regions 
(the blue areas in^a)). 

In conclusion, we have argued that the local gap in- 
homogeneity can be modeled by the modulation in the 
local superexchange interaction, induced by the dopant 
atoms. In the proposed strongly correlated electronic 
model, the local pairing interaction is controled by the 
superexchange interaction. Therefore, in the presence 
of the interstitial dopant atoms, the strongest effect of 
the doping would be an off-diagonal disorder that modu- 
lates superexchange on the neighboring bonds. Changes 
in the local superexchange lead to the local changes in 
the pairing strength and, in turn, in the local magnitude 
of the d-wave gap. An negative correlation of the hole 
density variation with the dopant location is also found. 
The self-consistent calculations also shows that the re- 
gions of large gap are registered to the dopant atoms. 
Large gap regions exhibit smaller and broader coherence 
peaks. These results are qualitatively consistent with re- 
cent STM observations by , such as positive corre- 
lation between large gap regions and oxygen dopant posi- 
tions. For a reasonable choice of parameters, we find the 
gap variations by about 20%. Although we use a specific 
RVB-like model, we believe the effect of local gap increase 
to be a general property of any model, where the pairing 
interaction is set by local superexchange. The following 
remarks are in order: (i) Our slave-boson mean-field cal- 
culation shows an increase of the average gap with a de- 
creased doping, where the dopant concentration becomes 
more dilute. The result is reasonably consistent with the 



experimental data [l9j. However, at this moment, we are 
unable to clarify whether the large gap in the low doping 
regime is dominated by the proximity to competing or- 
derings; (ii) Although the unscreened electrostatic poten- 
tials generated by the oxygen dopants do not tune the lo- 
cal pairing interaction in the proposed model, they might 
still play a role in the local electronic structure; (iii) To 
fully address various aspects of nanoscale electronic inho- 
mogeneity in high-Tc cuprates, a minimal effective model 
should take into account both the site-diagonal and off- 
diagonal disorder generated by the oxygen dopants, the 
resultant modulation in the superexchange interaction 
associated with the off-diagonal disorder, and possibly 
the residual Hubbard interaction. 
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